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Abstract. We perform an analytical study of the Integrated Sachs- Wolfe (ISW) 
effect within the framework of Unified Dark Matter models based on a scalar field 
which aim at a unified description of dark energy and dark matter. Computing the 
temperature power spectrum of the Cosmic Microwave Background anisotropics we are 
able to isolate those contributions that can potentially lead to strong deviations from 
the usual ISW effect occurring in a ACDM universe. This helps to highlight the crucial 
role played by the sound speed in the unified dark matter models. Our treatment is 
completely general in that all the results depend only on the speed of sound of the 
dark component and thus it can be applied to a variety of unified models, including 
those which are not described by a scalar field but relies on a single dark fluid. 
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1. Introduction 

Observations of large scale structure, search for la supernovae (SNIa), measurements 
of the Cosmic Microwave Background (CMB) anisotropics all suggest that two dark 
components govern the dynamics of the universe. In particular they are the dark matter 
(DM), responsible for structure formation, and an additional dark energy component 
that drives the cosmic acceleration observed at present [21 E]- Two main routes have been 
followed to provide a plausible realization of the dark energy, a non-zero cosmological 
constant A (see, e.g. Ref. [1]), or a dynamical dark energy (DE) component in the form 
of a scalar field, like for example Quintessence (U El El [TJ El EJ [101 EH] and k-essence 
models [T2"l [T3l [T3] . The latter are characterized by a Lagrangian with non-canonical 
kinetic term and inspired by earlier studies of k-inflation [15] [16] (a complete list of 
dark energy models can be found in the recent review [UJ). At different levels all of 
these scenarios suffer of non trivial fine tuning problems and of the so called cosmic 
coincidence problem (why Qdm and fl\ are both of order unity today). 

More recently the alternative hypothesis of unified models of dark energy and dark 
matter has been considered. In this case a single fluid behaves both as dark matter 
and dark energy. This has been variously referred to as "Unified Dark Matter" (UDM), 
or "Quartessence" . Among several models of k-essence considered in the literature 
there exist several UDM models. First of all purely kinetic models, i.e. the generalized 
Chaplygin gas (GCG) [HI [191 [20] model, the Scherrer and generalized Scherrer solutions 
[2T] |22j, a single dark perfect fluid with a simple 2-parameter barotropic equation of 
state [23], or the homogeneous scalar field deduced from the galactic halo space-times 
[21] . Moreover there are models in which imposing the Lagrangian of the scalar field to 
be a constant allows directly to describe a unified dark matter/ dark energy fluid [22] 
[25] (see also, with different approach, [26]). Alternative approaches to the unification 
of DM and DE have been proposed in Ref. [27] , in the frame of supersymmetry, and in 
Ref. [28] , in connection with the solution of the strong CP problem. 

The recent attitude in analysing the observational consequences of the DE models 
has been that of considering not only the background equation of state and its evolution 
with time, but also to focus on the sound speed which regulates the growth of the dark 
energy fluid perturbations on different cosmological scales. In this case the sound speed 
has been often treated as a completely independent parameter in order to explore the 
consequences on the CMB anisotropies and its effects on the low £ multipoles [29j EQl El] . 
The efficiency of this method relies on the observation that, for a single scalar field with 
canonical kinetic term, the speed of sound is equal to the speed of light, and thus it can 
cluster only on scales of the horizon size, while for other models it can be lower than 
unity, implying the possibility of clustering on smaller scales [32]. Another important 
issue is whether the dark matter clustering is influenced by the dark energy and, for 
the unified models, it becomes especially relevant in view of this approach. In the GCG 
model (both as dark energy and unified dark matter) strong constraints come from the 
CMB anisotropies [33l EH E2] and the analysis of the mass power spectrum |36j. In 
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the Scherrer solution the parameters of the model have to be fine-tuned in order for the 
model not to exhibit finite pressure effects in the non-linear stages of structure formation 

In this paper we consider cosmological models where dark matter and dark energy 
are manifestations of a single scalar field, and we focus on the contribution to the 
large-scale CMB anisotropics which is due to the evolution in time of the gravitational 
potential from the epoch of last scattering up now, the so called late Integrated Sachs- 
Wolfe (ISW) effect [38]. Through an analytical approach we point out the crucial 
role of the speed of sound in the unified dark matter models in determining strong 
deviations from the usual standard ISW occurring in the ACDM models. Our treatment 
is completely general in that all the results depend only on the speed of sound of the 
dark component and thus it can be applied to a variety of models, including those which 
are not described by a scalar field but relies on a single perfect dark fluid. In the case 
of ACDM models the ISW is dictated by the background evolution, which causes the 
late time decay of the gravitational potential when the cosmological constant starts to 
dominate [39]. In the case of the unified models there are two simple but important 
aspects: first, the fluid which triggers the accelerated expansion at late times is also 
the one which has to cluster in order to produce the structures we see today. Second, 
from the last scattering to the present epoch, the energy density of the universe is 
dominated by a single dark fluid, and therefore the gravitational potential evolution is 
determined by the background and perturbation evolution of just such a fluid. As a 
result the general trend is that the possible appearance of a sound speed significantly 
different from zero at late times corresponds to the appearance of a Jeans length (or a 
sound horizon) under which the dark fluid does not cluster any more, causing a strong 
evolution in time of the gravitational potential (which starts to oscillate and decay) 
and thus a strong ISW effect. Our results show explicitly that the CMB temperature 
power spectrum Ce for the ISW effect contains some terms depending on the speed of 
sound which give a high contribution along a wide range of multipoles i. As the most 
straightforward way to avoid these critical terms one can require the sound speed to be 
always very close to zero (thou see Sec. 3.2.3 for a more detailed discussion on this point). 
Moreover we find that such strong imprints from the ISW effect comes primarily from 
the evolution of the dark component perturbations, rather than from the background 
expansion history. The paper is organized as follows. In Sec. 2 we obtain the evolution 
equation for the gravitational potential. In Sec. 3 we start the analytical analysis of 
the ISW effect, dividing the resulting expression for the angular CMB power spectrum 
according to three relevant regions: those perturbation modes that enter the horizon 
after the acceleration of the universe becomes relevant, and perturbation modes that 
are inside or outside the sound horizon of the dark fluid. In Sec. 3.2.2 we point out 
those contributions to the ISW effect that are triggered by the sound speed and that are 
responsible for a strong ISW imprint. Sec. 4 contains our conclusions and a discussion 
of our results applied to various unified dark matter models. 



ISW effect in Unified Dark Matter Scalar Field Cosmologies: an analytical approach 4 
2. Linear perturbations in scalar field unified dark matter models 



(1) 



We consider the action that describes most of the dark matter unified models within 
the framework of k-essence 

S = Sq + Sip = J d^x^f^g 

where 

X = -^V^W • (2) 

We use units such that 8nG = 1 and signature (— , +, +, +). 
The energy-momentum tensor of the scalar field ip is 

^ = -;?=|^ = ^tr^^ + ^ *)*•■ < 3 > 

If X is time-like describes a perfect fluid with = (p+p)u^u u +p g^ v , with pressure 
C = p(<p,X), (4) 
and energy density 

p = p(p, X) = 2X^^1 - p(p, X) (5) 

where 

«, = ^ ■ (6) 

the four-velocity. 

Now we assume a flat, homogeneous Friedmann- Robertson- Walker background 
metric i.e. 

ds 2 = -dt 2 + a(t) 2 5 ij dx i dx j = a(r]) 2 (-dr] 2 + 5 ij dx i dx j ) , (7) 

where a(t) is the scale factor, o^- denotes the unit tensor and dr] — dt/a is the conformal 
time. In such the background evolution of the universe is characterized 

completely by the following equations 

U 2 = a 2 H 2 = l -a 2 p, (8) 

H'-H 2 = a 2 H = ~a 2 (p + p), (9) 

where Ti = a' /a, the dot denotes differentiation w.r.t. the cosmic time t and a prime 
w.r.t. the conformal time rj. On the background X = \p 2 = p' 2 /(2a 2 ) and the equation 
of motion for the homogeneous mode <p(t) becomes 

One of the relevant quantities for the dark energy issue is the equation of state w = p/p 
which in our case reads 

w = £ . (11) 
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On the other hand we will focus on the other relevant physical quantity, the speed 
of sound, which enters in governing the evolution of the scalar field perturbations. 
Considering small inhomogeneities of the scalar field 

<p(t,x)=po(t)+S<p(t,x), (12) 

we can write the metric in the longitudinal gauge as 

ds 2 = -(1 + 2<$>)dt 2 + (1 - 2<$>)a(t) 2 5 ij dx i dx j (13) 

since 5T t j = for i ^ j pED]. 

From the linearized (0 — 0) and (0 — i) Einstein equation one obtains (see Ref. [16] 



and Ref. [H]) 

V 2 $ = 1 - ^ ' rj [n^r + $) , (14) 



_. 2ih = l a 2 (p + p) ( kp_ 

2 ciu i/y 



2 $y la 2 (p + p) f^Scp 



a — 



^-7 + «f , (15) 



U) 2 H 2 v Vo 

where one defines a "speed of sound" c 2 relative to the pressure and energy density 
fluctuation of the kinetic term [TB] as 

_ (dp/dX) _ f 
s ~ (dp/dX) H + 2Xg| ' 1 ' 

Eqs. (1141) and f|T5l) are sufficient to determine the gravitational potential $ and the 
perturbation of the scalar field. Defining two new variables 

,is2 (^- ^KiH- <i7) 

where z = a 2 (p + p) 1 ^ 2 / (c s TC), we can recast and (TTo^) in terms of u and f jH] 

c a Au = zQV, CsU = 0Qy (18) 

where 9 = l/(c s z) = (1 + p/ p)- 1 ' 2 / (V?,a). Starting from jTB]) we arrive at the following 
second order differential equations for u [4T] 

m" - c^V 2 m - — u = . (19) 

Notice that this equation can also be used to describe any perfect fluid with equation of 
state p = p(p), up to a redefinition of c s . In this case c 2 = p'/ p' corresponds to the usual 
adiabatic sound speed. In this way, with the same equation (fT9~l) . we can also describe 
the ACDM model. Also pure kinetic Lagrangian (j4J) C(X) models (see for example 
Ref. |22j). can be described as a perfect fluid with the pressure p uniquely determined 
by the energy density, since they both depend on a single degree of freedom, the kinetc 
term X. 

Unfortunately we do not know the exact solution for a generic Lagrangian. However 
we can consider the asymptotic solutions i.e. long-wavelength and short- wavelength 
perturbations, depending whether c 2 k 2 -C \9"/9\ or c 2 k 2 ^> \9"/9\, respectively. This 
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means to consider perturbations on scale much larger or much smaller than the effective 
Jeans length for the gravitational potential Aj = c 2 \9j9"\. 

For a plane wave perturbation u oc Uk{i]) exp(ikx) in the short-wavelength limit 
(c 2 s k 2 > \6"/6\) we obtain 



CM 
cl /2 (r]) 



■ cos 



c s dr] 



(20) 



where C k is a constant of integration. Instead, neglecting the decaying mode, the long- 
wavelength solution (c 2 k 2 <C \9"/9\) is 

u k = A k (n)9 f_ , (21) 
where A k is a constant of integration . 

Once u is computed we can obtain the value of the gravitational potential $ through 
Eq. (|T7|) and the perturbation of the scalar field from Eq. (fi~5l) 



<fy> = 2V2X- 



a(p + p) 



(22) 



3. Analytical approach to the ISW effect 



Let us now focus on the ISW effect. The ISW contribution to the CMB power spectrum 
is given by 

2 



21 + 1 



ISW 



' " dk k 3 

4tt ~ 1 2tt 2 Jo T 21 + 1 

where 0} sw is the fractional temperature perturbation due to ISW effect 

k) = 2 jT tfft " « , (24) 

with i]q and 77* the present and the last scattering conformal times respectively and ji 
are the spherical Bessel functions. 

We now evaluate analytically the power spectrum (j23l . As a first step, following the 
same procedure of Ref. [39], we notice that, when the acceleration of the universe begins 
to be important, the expansion time scale 771/2 = rj(w = —1/2) sets a critical wavelength 
corresponding to krjin = 1. It is easy to see that if we consider the ACDM model then 
771/2 = ?7a i-e. when a\/ao = (^o/^a) 1 ^ 3 [39J. Thus at this critical point we can break 
the integral (!23l) in two parts [39J 



@isw 



(7/0, k) 



(23) 



21 + 1 

47T 



2tt 2 



Ietikfji/2 < 1) + ie,(^i/2 > 1) 



where 



and 



< 1) = / 

« 



1 /^i/ 2 rffc 
IT* 

K 



@ISW 



(??0, fc) 



21 + 1 



Ie t ( k Vi/2 > !) = / 



— fc 3 

l/»?l/2 K 



ep w (r/ ,A;) 



2/ + 1 



(25) 



(26) 



(27) 
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As explained in Ref. the ISW integrals (1241) takes on different forms in these two 
regimes 

Qzisw(7/o,fc) = f 2A$ fc ji[k(r] - r] 1/2 )} kr) 1/2 < 1 . . 

21 + 1 " \ 2& k ( Vk )I l /k kqi/2 > i 1 } 

where is the change in the potential from the matter-dominated (for example at 
recombination) to the present epoch 770 and rjk — t]q — (I + 1/2) / k is the conformal time 
when a given k-mode contributes maximally to the angle that this scale subtends on the 
sky, obtained at the peak of the Bessel fucntion j£. The first limit in Eq. ( [28]) is obtained 
by approximating the Bessel function as a constant evaluated at the critical epoch 7/1/2. 
Since it comes from perturbations of wavelenghts longer than the distance a photon can 
travel during the the time 771/2, a kick (2A§k) to the photons is the main result, and it 
will corresponds to very low multipoles, since 771/2 is very close to the present epoch 770 . 
It thus appears similar to a Sachs- wolfe effect (or also to the early ISW contribution). 
The second limit in Eq. (128]) is achieved by considering the strong oscillations of the 
Bessel functions in this regime, and thus evaluating the time derivative of the potentials 
out of the integral at the peak of the Bessel function, leaving the integral [39] 



With this procedure, replacing (f28h) in ( [26]) and (128b) in ( [271) we can obtain the ISW 
contribution to the CMB anisotropies power spectrum (!23l . 

Now we have to calculate, through Eqs. ([20l) - ([211 and ( TT7I) . the value of Q(k, 77) for 
kVi/2 ^ 1 and krji/ 2 3> 1. As we will see tha main differences (and the main difficulties) 
of the unified dark matter models with respect to the ACDM case will appear from the 
second regime of Eq. ( 1281) . 



3.1. Derivation of lQ t for modes kr\\j 2 < 1 

In the UDM models when ^771/2 <C 1 then c 2 s k 2 <C \9" ' /9\ is always satisfied. This is due 
to the fact that before the dark fluid start to be relevant as a cosmo logical constant, for 
V < V1/2, its sound speed generically is very close to zero in order to guarantee enough 
structure formation, and moreover the limit kr)i/ 2 1 involves very large scales (since 
77!/ 2 is very close to the present epoch). For the standard ACDM model the condition 
is clearly satisifed. In this situation we can use the relation ([21]) and becomes 



We immediately see that = $fc(0), the large scale gravitational potential during the 
radiation dominated epoch. The integral in Eq. ( 130]) may be written as follows 

a (r))dr] = I R + / a (fjjdfj , (31) 

•>r\R 
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where Ir = f^ R a 2 {rf)dfi and T}r is the conformal time at recombination. When 
r/i < 77 < tjr the UDM Models behave as dark matter Q. In this temporal range the 
universe is dominated by a mixture of "matter" and radiation and Ir becomes 



(r = V*a eq ( f + & + (32) 



where a eq is the value of the scalar factor at matter-radiation equality, £ = 77/77* and 
V* — (PeqO'l q /24)~ 1 / 2 = r] eq /{\/2— 1). With these definitions it is easy to see that an = 
a eq(&t + 2£r)- Notice that Eq. (1301) is obtained in the case of adiabatic perturbations. 
Since we are dealing with unified dark matter models based on a scalar field, there will 
always be an intrinsic non-adiabatic pressure (or entropic) pertubation. However for the 
very long wavelenghts, fc?7i/2 «C 1 under consideration here such an intrinsic perturbation 
turns out to be negligible [16J. For adiabatic perturbations ^k(vn) — (9/10)$fc(0) [40] 
and accounting for the primordial power spectrum, fc 3 |$fc(0)| 2 = Bk n ~ l , where n is the 
scalar spectral index, we get from Eq. (T28k ) 

IeAkVi/2 < 1) « 4(2/ + 1)73 J -k^jf^o - Vl/2 )] 



x 



1 H{ Vo ) 



a (fj)dfj 



(33) 



10 a 2 ( Vo ) 

where we have neglected Ir since it gives a negligible contribution. 

A first comment is in order here. There is a vast class of unified dark matter 
models that are able to reproduce exactly the same background expansion history of 
the universe as the ACDM model (at least from the recombination epoch onwards). For 
example this is the case of the the Scherrer and generalized Scherrer unified models [21] 
[22] , the generalized Chaplygin gas [HI [T9J, [20] for the parameter a which tends to zero, 
the models proposed in Ref. [22] and [25] where one impose the langrangian (i.e. the 
pressure) to be a constant, and also the model of a single dark perfect fluid proposed in 
Ref. [23]. For such cases it is clear that the low I contribution (|33|) to the ISW effect will 
be the same that is predicted by the ACDM model. This is easily explained considering 
that for such long wevelenght perturbations the sound speed in fact plays no role at all. 



3.2. Derivation of Iq 1 for modes krji/2 > 1 

As we have already mentioned in the previous section, in general a viable UDM must 
have a sound speed very close to zero for 77 < ^1/2 in order to behave as dark matter also 
at the perturbed level to form the structures we see today, and thus the gravitational 
potential will start to change in time for r\ > 771/2- Therefore for the modes £771/2 > 1, in 
order to evaluate Eq. (128b ) into Eq. fT2~T|) we can impose that r/^ > 771/2 which, from the 
definition of 77^ ~ 770 — (l+l/2)/k, moves the lower limit of Eq. ( |27l) to (1+ 1/2) /(rj — 7/1/2). 

I In fact the Scherrer [2T and generalized Scherrer solutions |22j in the very early universe, much 
before the equality epoch, have c s ^ and w > 0. However at these times the dark fluid contribution 
is sub-dominant with respect to the radiation energy density and thus ther is no substantial effect on 
the following equations. 
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Moreover we have that 771/2 
at the value of Defining 
V 



r] . We can use this property to estimate any observable 



X 



and 



K 



Vl/2 



kvi/2 , (34) 



we have 



taking a 



ak = a(Vk) = a(Xk) 



1, and 

dx 



(Xfc) = Vi/2&(Vk 



da 
dx 



*Xfc = 1 - Vi/2n , (35) 



Xo 



^x 2 



'/ + 1/2 N 



(36) 

Xo /v - XO ^ ' 

where 5\k = Xk — Xo — (Vk ~ Vo)/Vi/2 = —{I + 1/2)/ k. Notice that the expansion (1361) 
is fully justified, since as already mentioned above, the mimimum value of k in Eq. (1271) 
moves to (/ + 1/2)/ (770/771/2 ~~ 1)? making much less than 1. Therefore we can write 



I Qz isw (770, jOT 
(2Z + l) 2 



4 



4/f 



A- 



4/f 



(Xa 



c?X 



(Xo) 



dx 
d§ k 

dx d\ 2 



d 2 &k 

(xo)-rx(xo) 



'/ + 1/2 S 



+ 



K 



+ 



d 2 ^ 



d X 2 



(Xo) 



'1 + 1/2 S 



K. 



(37) 



In this case, during 771/2 < 77 < r] Q , there will be perturbation modes whose wave- 
length stays bigger than the Jeans length or smaller than it, i.e. we have to consider 
both the possibilities c 2 k 2 <C \0" /6\ and c 2 k 2 3> \9"/9\. In general the sound speed can 
vary with time, and in particular it might become significantly different from zero at 
late times. However, just as a first approximation, we exclude the intermediate situa- 
tion because usually 771/2 is very close to 770 (this situation will be briefly analyzed later). 



3.2.1. Perturbation modes on scales bigger than the Jeans length. 
When c 2 s k 2 <C \9"/9\ the value of $'(%) can be written from Eq.(j3"0]) as 

d 2 n /■* 



df^aL a(i)di ) 



(3? 



- 1 t=t(Vk) 

Now, using this expression in Eq. fl37j) . with the primordial power spectrum /c 3 |$fc(0)| 2 
_B/c n_1 , the value of ff2~7]) may be written as 

2 



WM1/2 > 1 
2/ + 1 



4/^</- 2 1 



1+1/2 K 3 



K 



n-l 



4I 2 Brt 



/2 



1 / Xo - 1 
3-77 W + 1/2 



dx 

\ 3— n 



(Xk) 



d<h 

dx 



(Xo) 



+ 
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2(2 + 1/2) fxo-1 



+ 



4-71 

(/ + 1/2) 2 



n 



v Z + l/2, 

' xo-l N 
1 + 1/2, 



■(xo)-r-r(xo)+ 



2 



(Xo) 



(39) 



with (d$fc/^x)xo = ^i/2$' fc (^o) and with (d 2 ^ k /dx 2 



>X(> 



A second relevant comment follows from the fact that If ~ 1/2 for I 3> 1. We thus see 
that for n — 1 and for / ^> 1 the contribution to the angular power spectrum from the 
modes under consideration is l(l + l)C( sw /(4tt) = Z(Z + l)I e ,(A»7i/2 > l)/(2vr 2 (2/ + l)) ~ 



1/Z. In other words we find a similar slope as in [391 H2] found in the ACDM model. 
Recalling the results of the previous section, this means that in the unified dark matter 
models the contribution to the ISW effect from those perturbations that are outside the 
Jeans length is very similar to the one produced in a ACDM model. The main difference 
on these scales will be present if the background evolution is different from the one in 
the ACDM model, but for the models where the background evolution is the same, as 
those proposed in Refs. [211 E2J, [251 E3] no difference at all can be observable. 



3.2.2. Perturbation modes on scales smaller than the Jeans length. 

When c 2 s k 2 3> \Q" /0\ we must use the solution (TSUI) and through the relation fflTk) the 
gravitational potential is given by 

ft] 



Mv) = \[(p + p)/c s } 1/2 (v)c k (m 



12) cos I k I c s (fj)dfj . (40) 

V •''71/2 / 



In Eq. (14"0]) Cfc (771/2) = < ^fc(0)Ci/ 2 is a constant of integration where 



Ci 



/2 



(41) 

[(p + p)/c s ] 1/2 (771/2) 

and it is obtained under the approximation that for 77 < 771/2 one can use the 
longwavelenght solution ( l30l) . since for these epochs the sound speed must be very close 
to zero. Notice that Eq. ([40]) shows clearly that the gravitational potential is oscillating 
and decaying in time. 

Defining for simplicity C = C 2 , 2 [(p + p)/c s ] (770) /4, we take the time derivative 
of the gravitational potential appearing in Eq. (128b ) by employing the expansion of 
Eq.CT. We thus find that Eq. (ETJ) yields 
Ie,(kr) 1/2 > 1) 



21 + 1 



AC* Blfa^ 1 \ C {k5>l2>c 2 } (l + 1/2) 
00 dn 



dn 



+ 



+ C{k4,ii, c 2}(l + I/ 2 ) 
(Z + l/2) 2 



Jl+Hl K V U ' 



+ 



00 ck 



1+1/2 
XQ-1 



k" 1 cos ( Do k) sin(D K) 



+ 



+ 



C{k3,l0,c 2 } + C{fc3,«2,c 2 }(2 + 1/2)" 



dn 



.n-l 



1+1/2 
xo-i 



COS (DnK,) 



K 



+ 
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f(P + P),xx f(P + P),x 
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C{k2,l0,sc} 

C 



-4 



{fcUO.s 2 } — 4c s | XQ , L-{fcl,Z2,c 2 } — ^ c sl X0 ) ^{fel,U,ac} — si x „ 

In this case we have defined (-) >x = d(-)/d\ and we recall that the dimensionless 
variables x an d k are defined in Eq. (|34l) . We have indicated the coefficients 
^{fcbVMiM) in such a way to signal that they multiply an integral in k of K n ~ l j k? 
times sin(Do K ) cos(_Do K ) and the overall multipole coefficent is (I + l/2)\ We can infer 
from (j4"2"j) that for n < 1 all integrals are convergent. Notice that a natural cut-off 
in the various integrals is introduced for those modes that enter the horizon during 
the radiation dominated epoch, due to the Meszaros effect that the matter fluctuations 
will suffer until the full matter domination epoch. Such a cut-off will show up in the 
gravitational potential and in the various integrals of Eq. f|42l) as a (k eq /k) A factor, where 
k eq is the wavenumber of the Hubble radius at the equality epoch. 

A simple inspection of Eq. (142]) shows one of our main results. The terms of Eq. (l4"2{) 
where the coefficients C turn out to be proportional to the sound speed c s cause the 
growth of /(/ + l)/ei(^ ? 7i/2 > 1)/(2Z + 1), (and hence of the power spectrum 1(1 + 1)C; 
through Eq. (l25l) ). as / increases. This means that, if the sound speed of the unified 
dark matter fluid starts to differ significantly from zero at late times, the consequence 
is to produce a very strong ISW effect, and clearly this does not happen in a ACDM 
universe since c 2 s = always. This effect is easily explained by considering that the 
energy density of the universe in the unified models is dominated at late time by a just 
single fluid. Therefore an eventual appearance of a Jeans lenght (i.e. a departure of the 
sound speed from zero) makes the oscillating behaviour of the dark fluid pertubations 
under the Jeans lenght immediately visible through a strong time dependence of the 
gravitational potential. In fact one can verify that the scalar field fluctuations (|22l) are 

1 /2 

oscillating and decaying in time as 5(p ~ (k/a) [c s / (dp/dX)] ' sin(fc ffj c s dn). Similar 
results have been discussed in the case of the GCG model in Refs. [3H 135] . 

We point out that the potentially most dangerous term in Eq. (1421) is the one 
identified by the coefficient C^i,i2,c 2 } 



(P + P),x 
(P + P) 
C ( 
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Xp + p),x 
ip + p) 



8ci 



(43) 



1 + 1/2 

xo-i 



-K n 1 COS' 



(44) 



Such a term makes the power spectrum l(l + l)Ci to scale as Z 3 until / ~ 25. This angular 
scale is obtained by considering the peak of the Bessel functions in correspondence of 
the cut-off scale k eq , I rs k eq (r] — 771/2) ■ In fact, for smaller scales, the integral identified 
by the coefficient C{ fcl / 2 ,c 2 } will decrease as l/£. 



3.2.3. Intermediate case. 

Now we shall briefly discuss the intermediate case that corresponds to perturbation 
modes that initially are outside the Jeans length and then, due to a time variation 
of the sound speed, they fall inside them. This corresponds to consider the range 
[(/ + l/2)/(xo - I)] 2 < «?j = \B, XX /B\ /(* < kL. In this case k > (1+ l/2)/( Xo - 1) and so 
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we can use the same procedure described before. Indeed when k ~ kj i.e. c 2 s k 2 j ~ \0"/8\ 
it can be written as follows 

/Z + l/2\l 2 1 



c 2 s \o, xx /or - c 2 s \e, xx /9\- : > 

XO 



XO 



(45) 



i.e. 



K j = Kj (l) = B t + yjBf + A (46) 

where B x = {[ln(c^)] )X - [In \9, xx /9\\, x }\ x pl + l)/4 and A = [|0, xx /0| /(2^)]| xo . We 
immediately see that IeXkVi/2 > 1) can be divided into two parts. The first part is 
identical to fl39|) except for the upper limit of the integral. Indeed now the upper limit 
is Kj(l). In order to derive the second part we note that now the lower limit of the 
integral in k is kj(1) and that Uk(r]k) = [Ck(i]j)/ c l^(Vk)] cos [k J^* c s (fj)dffj where r/j is 
the conformal time when c 2 \ rij k 2 ~ \6"/6\ VJ (Vj * s function of k) and 

CM = 2 ** (0) 1 ZtvJrJ ' (47) 

4. Discussion of some examples and conclusions 

In most of the UDM models there are several properties in common. It is easy to see 
that in Eq. (j3"Tl) Ir is negligible because of the low value of a eq . 

Moreover in the various models usually we have that strong differences with respect to 
the ISW effect in the ACDM case can be produced from those scales that are inside the 
Jeans length as the photons pass through them. For these scales (which depend on the 
particular model) the perturbations of the unified dark matter fluid play the main role. 
On larger scales instead we find that they play no role and ISW signatures different 
from the ACDM case can come only from the different background expansion history. 
We have found that when k 2 ^> k 2 = c~ 2 \9" /9\ (see (fT9j) ) one must take care of the 
term in Eq. (f4"2|) proportional to C^i,i2,c 2 }- Indeed this term grows faster than the other 
integrals contained in fj4~2]) when / increases up to / ~ 25. It is responsible for a strong 
ISW effect and hence it will cause in the CMB power spectrum /(/ + 1)C// (2ti) a decrease 
in the peak to plateau ratio (once the CMB power spectrum is normalized). In order to 
avoid this effect, a sufficient (but not necessary) condition is that all the models have 
to satisfy c 2 k 2 < \9"/9\ for the scales of interest. The maximum constraint is found in 
correspondence of the scale at which the contribution Eq. (j4~2l proportional to C{ fcl : i2,c 2 } 
takes it maximum value, that is k w k eq . For example in the Generalized Chaplygin Gas 
model (GCG), i.e when p = — A 1 ^ 1 ^ / p a and c 2 = —aw, we deduce that \a\ < 10~ 4 
(see Refs. j2U] [M] [22] and [13]). This is also in accordance with [36] which performs an 
analysis on the mass power spectrum and gravitational lensing constraints thus finding 
a more stringnet constraint. 

As far as the generalized Scherrer solution models [22] are concerned, in these models 
the pressure of the unified dark matter fluid is given by p = g n (X — X ) n — A, where g n 
is a suitable constant and n > 1. The case n = 2 corresponds to unified model proposed 
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by Scherrer [2T]. In this case we find that imposing the constraint c 2 k 2 < \9"/9\ for the 
scales of interest we get that e = (X — X Q )/X < (n — 1) 1CT 4 . 

If we want to study in greater detail what happens in the GCG model when c 2 s k 2 ^> \9" /9\ 
we discover the following things: 

• for 10 -4 < a < 5 x 1CT 3 , where we are in the "Intermediate case". Now c 2 = —aw 
is very small and the background of the cosmic expansion history of the universe 
is very similar to the ACDM model. In this situation the pathologies, described 
before, are completely negligible. 

• When treating 6 x 10~ 3 < a < 1 a very strong ISW effect will be produced and 
we have estimated the same orders of magnitude for the decrease of the peak to 
plateau ratio in the anisotropy spectrum 1(1 + 1)C[/(2tt) (once it is normalized) 
that can be inferred from the authors of [M] obtained in the numerical simulations 
(having assumed that the production of the peaks during the acosutic oscillations 
at recombination is similar to what happens in a ACDM model, since at recombi- 
nation the effects of the sound speed will be negligible). 

An important observation arises when considering those UDM models that reproduce 
the same cosmic expansion history of the universe as the ACDM model. Among these 
models one can impose the condition w = —c 2 which, for example, is predicted by UDM 
models with a kinetic term of the Born-Infeld type [22] [25] [26]. In this case, computing 
the integral in Eq. (f4"2]) proportional to C{fci,z2,c 2 } which give the main contribution to 
the ISW effect we have estimated that the corresponding decrease of peak to plateau 
ratio is about one third with respect to what we have in the GCG when the value of a 
is equal to 1. The special case a = 1 is called "Chaplygin Gas" (see for example [19] ) 
and it is characterized by a background equation of state w which evolves in a different 
way to the standard ACDM case. 

From these considerations we deduce that this specific effect stems only in part from the 
background of the cosmic expansion history of the universe and that the most relevant 
contribution to the ISW effect is due to the value of the speed of sound c 2 . 
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